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Abstract. Relic density calculations are often used to constrain particle physics models, and in particular supersymmetry. 
We will show that the presence of additional energy or entropy before the Big-Bang nucleosynthesis can however completely 
change the relic density constraints on the SUSY parameter space. Therefore one should be extremely careful when using the 
relic density to constrain supersymmetry as it could give misleading results, especially if combined with the future collider 
data. Alternatively, we will also show that combining the discoveries of the future colliders with relic density calculations 
can shed light on the inaccessible pre-BBN dark time physics. Finally we will present Superlso Relic, a new relic density 
calculator code in Supersymmetry, which incorporates alternative cosmological models, and is publicly available. 
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INTRODUCTION 

The recent cosmological observations give evidence for the presence of a cosmological matter density, which repre- 
sents about 27% of the total density of the Universe, and the presence of the so-called dark energy. Using the total 
matter density observed by WMAP [1] and the baryon density indicated by Big-Bang nucleosynthesis (BBN) [2, 3] 
and including the theoretical uncertainties, the dark matter density range at 95% C.U. is deduced [4]: 

0.094 <£l DM h 2 < 0.135, (1) 

where h is the reduced Hubble constant. In the following, we also refer to the older range [5] 

0.1 < £l DM h 2 < 0.3 . (2) 

In supersymmetric models the lightest stable supersymmetric particle (USP) constitutes the favorite candidate for dark 
matter. If the relic density can be calculated precisely, the accuracy of the latest WMAP data can therefore be used 
to constrain the supersymmetric parameters. The computation of the relic density is well-known within the standard 
model of cosmology [6, 7], and is implemented in automatic codes, such as MicrOMEGAs [8], DarkSUSY [9] or 
Superlso Relic [10]. 

Nevertheless, the nature of the dark energy and the properties of the Universe in the pre-BBN epoch are still 
unknown, and the BBN era, at temperatures of about 1 MeV, is the oldest period in the cosmological evolution when 
reliable constraints are derived. The cosmology of the primordial Universe could therefore be much more complex, 
and the pre-BBN era could have for example experienced a slower or faster expansion. Such a modified expansion, 
even though still compatible with the BBN or the WMAP results, would change the LSP freeze-out time and the 
amount of relic density (see for example [4, 11, 12, 13]). A similar question exists concerning the entropy content of 
the Universe at that period, and the eventual consequences on the energy conservation in the primodial epochs (see for 
example [14, 15, 16, 17]). 

In the following, we will first describe the way the relic density is calculated in the standard cosmology. We will 
then present how the calculations can be modified in altered cosmological scenarios. We then analyze the consequences 
of the cosmological uncertainties on the supersymmetric parameter searches. We also inverse the problem and show 
that the determination of a beyond the Standard Model particle physics scenario will give hints on the cosmological 
properties of the early Universe. Finally, we will present the Superlso Relic package and conclude. 
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RELIC DENSITY IN STANDARD COSMOLOGY 



The cosmological standard model is based on a Friedmann-Lemaitre Universe, approximately flat, incorporating a 
cosmological constant accelerating its expansion, and filled with radiation, baryonic matter and cold dark matter. 
Before BBN, the Universe expansion is dominated by the radiation density, and therefore the expansion rate H of the 
Universe is determined by the Friedmann equation 

H 2 = *-f Prad , (3) 

where 

Pra d (T)=g ef f(T)^T 4 (4) 

is the radiation density and g Q ff is the effective number of degrees of freedom of radiation. The computation of the 
relic density is based on the solution of the Boltzmann evolution equation [6, 7] 

dn/dt = —3Hn — ((7 e ffv) (n 2 — ^q) > (5) 

where n is the number density of all supersymmetric particles, n&q their equilibrium density, and (<7 e ffv) is the thermal 
average of the annihilation rate of the supersymmetric particles to the Standard Model particles. By solving this 
equation, the density number of supersymmetric particles in the present Universe and consequently the relic density 
can be determined. 

The computation of the thermally averaged annihilation cross section (o" e ffv) requires the computation of the many 
annihilation and co-annihilation amplitudes. The annihilation rate of supersymmetric particles / and j into SM particles 
k and / is defined as [6, 7]: 

where jtft is the transition amplitude, s the center-of-mass energy, gt the number of degrees of freedom of the particle 
i and the final center-of-mass momentum such as 

[s - (m k + mi) 2 ] 1/2 [s - (m k - mi) 2 ] 1/2 

Pu = 7t-f • y') 

m Eq. (6) is a symmetry factor equal to 2 for identical final particles and to 1 otherwise, and the integration is over 
the outgoing directions of one of the final particles. Moreover, an average over initial internal degrees of freedom is 
performed. 

The effective annihilation rate W e & can be defined by 

glspPesWeK = L SigjPij^ij (8) 

ij 

with 



Pe&(Vs) = ^\/(V^) 2 -4m 
The following relation can therefore be deduced: 
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where is the angle between particles i and k. The thermal average of the effective cross section is then obtained by: 
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where K\ and K 2 are the modified Bessel functions of the second kind of order 1 and 2 respectively. 
The ratio of the number density to the radiation entropy density Y(T) = n(T)/s(T) is defined, in which 



IK 2 a 

s(T) = h eS (T) — T 3 . (12) 

h Q ff is the effective number of entropic degrees of freedom of radiation. Combining Eqs. (3) and (5) and defining the 
ratio of the LSP mass over temperature x = mj^sp / T, yield 



dY 



dx V 45G 

with 



1/2 

! M(^-Iq)> d3) 



*y 2 =^(i+^^i, a4) 



/Self V 3/ *eff Jr 
and 

>W = ^^E^2(y) . (15) 



where / runs over all supersymmetric particles of mass m ; and with g; degrees of freedom. The freeze-out temperature 
Tf is the temperature at which the LSP leaves the initial thermal equilibrium when Y(Tf) = (1 + S)Yeq(7/)» with 
5 ~ 1.5. The relic density is obtained by integrating Eq. (13) from x = to m^p/To, where To = 2.726 K is the 
temperature of the Universe today [6, 7]: 

n LSP /, 2 = ^§P«jW w 2 755 x 10^F(7b) , (16) 
where is the critical density of the Universe, such as 

Ht=*f P «, (17) 

Hq being the Hubble constant. 



RELIC DENSITY IN ALTERNATIVE COSMOLOGICAL SCENARIOS 

In presence of non-thermal production of SUSY particles, the Boltzmann equation becomes 

^ = -3Hn- {av)(n 2 -n 2 eq )+N D . (18) 

The term No is added to provide a parametrization of the non-thermal production of SUSY particles. The expansion 
rate H can also be modified: following [4, 17], po is introduced as an effective dark density which parametrizes the 
expansion rate modification. The Friedmann equation then becomes 

o 87TG 

H 2 = —(p md +p D ), (19) 

where p ra j is the radiation energy density, which is considered as dominant before BBN in the standard cosmological 
model. In case of additional entropy fluctuations, the entropy evolution reads 

^ = -3tfj + £z), (20) 
at 

where s is the total entropy density. T^d parametrizes here effective entropy fluctuations due to unknown properties of 
the early Universe. 
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Separating the radiation entropy density from the total entropy density, i.e. setting s = s rac i where s rac i is the 
radiation entropy density and sd is an effective entropy density, the following relation between sd and can be 
derived: 

1 h G f£ _ dso 



4tt 3 G 



(21) 



Following the standard relic density calculation method [6, 7], we introduce Y = n/s, and Eq. (18) becomes 



dY 

dx 
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where x = m L sp/T, m L sp being the mass of the relic particle, 



sd 



2k 1 



and 



A effC r ) 45 

45 



7 3 



Pd 



(22) 



7r£ T 4 

on J- 



The relic density can then be calculated in the standard way: 

Qh 1 = 2.755 x 10 8 7 ^zWGeV . 



(23) 



(24) 



(25) 



where Yo is the present value of Y. In the limit where p D = sd = ^d = Nd = 0, usual relations are retrieved. We should 
note here that sd and are not independent variables. In the following, we neglect 
We use the parametrizations described in [4, 17] for p D and sd- 



PD = KpPrad(T B BN) f J— J 



and 



5d = K s S r ad(TBBN) 



f t y s 

\Tbbn J 



(26) 



(27) 



where T BB n is the BBN temperature. K p (k s ) is the ratio of effective dark energy (entropy) density over radiation 
energy (entropy) density at BBN time, and n p and n s are parameters describing the behavior of the densities. We refer 
to [4, 17] for detailed descriptions and discussions on these parametrizations. 



SUPERSYMMETRIC CONSEQUENCES 

Relic density is often used to constrain SUSY parameter space (see for example [18]). In particular the non-universal 
Higgs model (NUHM) provides attractive candidates for dark matter [19]. For our analysis, we consider the NUHM 
parameter plane (jU,mA), fixing the other parameters (mo = 1 TeV, m ly / 2 = 500 GeV, tan/3 = 35, Ao = 0). About 
250,000 random SUSY points in the NUHM parameter plane (jijha) are generated using SOFTSUSY v2.0.18 [20], 
and for each point we compute flavor physics observables, direct limits and the relic density with Superlso Relic v2.7 
[10]. In Fig. 1, the excluded zones due to different observables are displayed: the red area is excluded by the isospin 
asymmetry of B — > K*y, the green by the inclusive branching ratio of b — > sy, the yellow area leads to tachyonic 
particles and the gray zone has been excluded by collider searches. All these exclusions are related to particle physics 
and are subject to uncertainties which are under control. The dark (light) blue zones are favored by the WMAP (old) 
dark matter constraints. Hence, in the top left figure for example, which corresponds to the standard cosmological 
model, only tiny strips remain compatible with all the displayed constraints, and the relic density observable happens 
to be extremely constraining. This kind of figures is often used to determine the favored SUSY parameter space zones. 
In the next plots of Fig. 1 however, we show that the presence of additional energy density in the early Universe, even 
negligible at BBN time, can completely change the results. These plots show the influence of a quintessence-like dark 
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FIGURE 1. Constraints in the NUHM parameter plane (/^m^) in presence of additional dark energy for several values of K p 
and n p , and for mo = 1 TeV, m^ = 500 GeV, tan (5 = 35, Aq = 0. The color code is given in the text. 

energy (n p = 6) whose proportions relative to the radiation density at BBN time are respectively K p = 10 -5 , 10 -4 , 
10 -3 and 10 -2 . The last plot shows the influence of a density of a decaying scalar field (n p = 8) with an extremely 
low K p = 10 -5 . In these plots the relic favored zones are moved towards the center, completely modifying the favored 
SUSY parameters. 

The result is similar when considering the influence of additional entropy. In Fig. 2, the first plot shows for reference 
the constraints in the standard model of cosmology. The second plot presents the influence of a dark entropy density 
with n s = 4 and K s = 10 -3 , which can occur in the case of reheating. The next plots show the influence of a dark 
entropy density with n s = 5 and k s — 10 -4 , K*^ = 10 -3 and k s — 10 -2 . The relic density favored areas are this time 
moved outwards with different shapes. 

From these two analyses, it is clear that the relic density calculations can be strongly altered by the presence of 
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FIGURE 2. Constraints in the NUHM parameter plane (/x, tha) in presence of additional dark energy for several values of K s and 
n s , and for mq = 1 TeV, mi/ 2 = 500 GeV, tan /3 = 35, Aq = 0. The color code is given in the text. 

entropy or energy densities, even very small or negligible at the BBN time (i.e. completely unobservable in the current 
cosmological data). Thus, cosmological unknown properties of the early Universe can completely change the favored 
SUSY parameter space, and could lead to erroneous affirmations on the properties of the SUSY particles. 



INVERSE PROBLEM 



We have seen in the previous section that it is critical to use the relic density to constrain SUSY. However the problem 
can be inverted: the future colliders will give the possibility to determine the new physics and relic particle properties 
[21], and by calculating the relic density in different cosmological scenarios it would be possible to determine or to 
verify some of the physical properties of the early Universe. In such cases, combining the relic density constraints with 
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FIGURE 3. Constraints in the NUHM parameter plane (fi, m&) in presence of additional dark entropy and energy for K p = 10 2 , 
n p = 6, K s = 10 -2 , n s = 5 (left) and K p = 10~ u , n p = 8, K s = 10 -2 , n s = 4 (right). The color code is given in the text. 
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FIGURE 4. Constraints from Y p (left) and 2 H/H (right) on the dark energy parameters (n p , Kp). The parameter regions excluded 
by BBN are located between the black lines for Y p , and in the upper left corner for 2 H/H. The colors correspond to different values 
of Y p and 2 H/H. 



the BBN limits would remove some of the degeneracies. In Fig. 3 for example, the constraints are presented for two 
different cosmological scenarios: K p — 10 -2 , n p = 6, K s = 10 -2 , n s = 5 on the left, K p = 10 -11 , n p = 8, K s = 10 -2 , 
n s = 4 on the right. The cosmological properties are different, but the obtained zones are rather similar, which means 
that there is a degeneracy between dark energy and dark entropy effects. Let us assume that the subjacent SUSY model 
leads to a relic density inside the WMAP favored zone of these plots, and disfavored in the standard cosmology model. 
In Figs. 4 and 5 the current limits on the dark energy and entropy properties from the Y p and 2 H/H BBN constraints 
are presented. The areas on the top of the black lines lead to unfavored element abundances. Comparing Fig. 3 with 
Figs. 4 and 5, we notice that the left plot (K p = 10 -2 , n p = 6, K s = 10 -2 , n s = 5) is disfavored by BBN, while the 
right one is still compatible. Thus, in such cases, knowing the particle physics properties will enable us to recognize 
whether the favored cosmological model is the standard one, and by considering also the BBN data, we can already 
distinguish between different cosmological scenarios. 



SUPERISO RELIC 

The results presented here are obtained using the public code Superlso Relic v2.7 [10], available on 
http://superiso.in2p3.fr/relic. Superlso Relic is an extension of the Superlso program [22], which is 
devoted to the calculation of flavor physics observables in the Two-Higgs Doublet Model, Minimal Supersymmetric 
Standard Model and Next-to-Minimal Supersymmetric Standard Model. The main purpose of Superlso Relic is to 
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FIGURE 5. Constraints from Y p (left) and 2 H/H (right) on the dark entropy parameters (n s , K s ). The parameter regions excluded 
by BBN are located in the upper left corner for Y p , and in both upper corners and in between the right black lines for 2 H/H. The 
colors correspond to different values of Y p and 2 H/H. 



compute the relic density in the cosmological standard model as well as in alternative scenarios. All models described 
here are already implemented, and other models will soon be included. For more information on Superlso Relic, we 
refer to its website, or to the manual [10]. 

CONCLUSION 

We showed that using the relic density to constrain the supersymmetric parameter space is highly dependent on the 
cosmological assumptions, and it is then critical to use the relic density constraints for the scientific preparation of the 
future colliders. We also noticed that inverting the problem, i.e. using the particle physics results to determine the early 
Universe properties, can be of great interest as it will give access to an unknown part of the Universe history. 
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